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SUMMARY 


Inverse  synthetic  aperture  radar  (ISAR)  is  an  imaging  technique  that  shows  great 
promise  in  classifying  airborne  targets  in  realtime  under  all  weather  conditions.  Over  the 
past  several  years,  a  large  body  of  ISAR  data  has  been  collected  and  considerable  effort 
has  been  expended  to  develop  algorithms  to  form  high-resolution  images  from  these 
data.  One  important  goal  of  workers  in  this  field  is  to  develop  software  that  will  do  the 
best  job  of  imaging  under  the  widest  range  of  conditions.  The  success  of  classifying 
targets  using  ISAR  is  predicated  upon  forming  highly  focused  radar  images  of  these 
targets.  Efforts  to  develop  highly  focused  imaging  computer  software  have  been 
challenging,  mainly  because  the  imaging  depends  on  and  is  affected  by  the  motion  of  the 
target,  which  in  general  is  not  known.  Specifically,  the  target  usually  has  both  rotational 
motion  about  some  axis  and  translational  motion  as  a  whole  with  respect  to  the  imaging 
radar.  The  slant-range  translational  motion  kinematic  quantities  of  the  target  must  be 
first  accurately  estimated  from  the  radar  data  and  compensated  before  a  focused  image 
can  be  formed.  Following  slant-range  translational  motion  compensation,  the  image  is 
further  focused  by  determining  and  correcting  for  target  rotation. 

We  propose  the  use  of  the  burst  derivative  measure  as  a  means  to  improve  the 
computational  efficiency  of  currently  used  motion  compensation  ISAR  algorithms.  This 
technical  report  describes  the  use  of  the  burst  derivative  measure  in  motion 
compensation  ISAR  algorithms  for  estimating  both  the  slant-range  translational  motion 
and  rotational  kinematic  quantities  of  an  uncooperative  airborne  target.  Tests  have  been 
performed  on  simulated  as  well  as  actual  ISAR  data  using  both  a  Sun  4  SPARCstation  2 
and  a  parallel  processing  transputer  array.  Preliminary  results  indicate  that  the  burst 
derivative  measure  gives  significant  improvement  in  processing  speed  over  the 
traditional  entropy  measure  presently  employed. 
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1.0  INTRODUCTION 


Potentially  the  most  powerful  means  to  rapidly  classify  airborne  targets  at  long  range 
under  all  weather  conditions  is  through  the  use  of  some  type  of  imaging  radar 
(Ausherman  et  al.,  1984).  Inverse  synthetic  aperture  radar  (ISAR)  is  one  technique  that 
shows  real  promise  in  this  area  (Wehner,  1987).  Scientists  and  engineers  at  the  RDT&E 
Division  of  the  Naval  Command,  Control  and  Ocean  Surveillance  Center  (NCCOSC)  in 
San  Diego  have  developed  microwave  stepped-frequency  ISAR  systems  as  well  as  the 
algorithms,  software,  and  processing  architectures  necessary  to  process  ISAR  signal 
returns  for  displaying  highly  focused  ISAR  images. 

The  translational  motion  compensation  (TMC)  algorithm,  one  of  two 
computationally  intensive  ISAR  algorithms  presently  employed,  estimates  the  slant- 
range  translational  velocity,  acceleration,  and  jerk  of  an  uncooperative  target  with 
respect  to  the  imaging  radar.  The  determination  of  the  target's  translational  kinematic 
information  is  the  first  important  step  toward  achieving  a  highly  focused  ISAR  image. 
The  input  to  the  TMC  algorithm  is  a  two-dimensional  unprocessed  ISAR  spatial 
frequency  spectrum.  This  spectrum  is  multiplied  by  a  spatial  frequency  filter  containing 
initial  guess  information  of  the  target's  slant-range  translational  velocity,  acceleration, 
and  jerk.  The  resulting  filtered  spectrum  is  then  Fourier  transformed  using  the  two- 
dimensional  fast  Fourier  transform  (2D-FFT)  algorithm.  The  sharpness  of  the  resulting 
space  domain  ISAR  image  is  determined  by  computing  its  entropy.  This  process  is 
repeated  a  number  of  times  using  different  values  of  velocity,  acceleration,  and  jerk  in 
the  spatial  frequency  filter.  Those  values  that  give  rise  to  an  ISAR  image  having 
minimum  entropy  are  used  as  the  best  estimates  of  the  target's  actual  slant-range 
translational  velocity,  acceleration,  and  jerk.  Once  the  translational  kinematic  estimates 
have  been  determined,  then  processing  may  advance  to  the  second  stage  where  the 
rotational  kinematics  of  the  target  are  estimated  using  the  rotational  motion 
compensation  (RMC)  algorithm.  The  RMC  algorithm  is  used  to  further  improve  the 
sharpness  of  the  translational  motion-compensated  ISAR  image  obtained  as  the  output 
of  the  TMC  algorithm. 

The  authors  have  been  investigating  the  use  of  a  new  measure,  which  we  call  the 
burst  derivative  measure,  as  an  alternative  means  for  estimating  the  slant-range 
translational  and  rotational  kinematic  quantities  of  an  airborne  target.  The  burst 
derivative  measure  is  an  outgrowth  of  an  earlier  NCCOSC  RDT&E  Division  attempt  in 
using  the  Fisher  measure  (Frieden,  1989, 1990)  as  a  means  to  improve  the  computational 
throughput  rates  of  the  TMC  and  RMC  algorithms. 
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The  burst  derivative  measure  offers  a  significant  speed  advantage  over  the  entropy 
measure  in  that  computations  may  be  confined  strictly  to  the  spatial  frequency  domain. 
Hence,  the  numerous  2D-FFTs  presently  required  using  the  entropy  measure  in  the  TMC 
and  the  RMC  algorithms  could  be  avoided.  This  technical  report  describes  preliminary 
results  on  the  use  of  the  burst  derivative  measure  in  estimating  the  translational  and 
rotational  kinematics  of  an  uncooperative  airborne  target. 

In  the  following  sections,  we  will  present  the  underlying  theory  pertinent  to  the 
formation  of  the  motion-compensated  ISAR  image.  We  first  begin  our  discussion  with  a 
description  of  the  ISAR  geometry,  the  formation  of  the  unprocessed  ISAR  spatial 
frequency  spectrum,  a  description  of  the  TMC  and  RMC  algorithms,  and  finally 
preliminary  results  for  both  simulated  as  well  as  real  aircraft  data. 


2.0  ISAR  GEOMETRY 


The  geometry  employed  for  the  ISAR  scenario  is  depicted  in  figure  1.  For 
pedagogical  purposes,  lets  consider  a  two-dimensional  aircraft  moving  in  a  plane.  A 
stationary  ground-based  radar  (assumed  to  be  a  monostatic  omnidirectional  transmitter 
and  receiver  for  simplicity)  is  located  at  the  origin  of  a  right-handed  Cartesian 
coordinate  system  XY.  A  solid  extended  aircraft  is  moving  relative  to  the  ground-based 
radar.  The  center  of  mass  of  this  aircraft  is  located  at  the  origin  of  a  second  right-handed 
Cartesian  coordinate  system  xy  whose  corresponding  axes  are  parallel  to  those  of  XY.  In 
addition,  a  third  right-handed  Cartesian  coordinate  system  x'y',  attached  to  the  aircraft, 
is  used  in  describing  the  rotation  of  the  aircraft  with  respect  to  the  xy  system.  The  origins 
of  xy  and  x'y'  are  coincident.  The  angular  displacement,  0  (t) ,  describes  the  rotational 
position  of  the  aircraft  relative  to  xy.  The  slant-range,  r(t),  describes  the  translational 
position  of  the  aircraft's  center  of  mass  relative  to  the  radar.  Recall  from  classical 
mechanics  that  the  overall  motion  of  an  extended  rigid  body  can  be  described  by  the 
translational  motion  of  its  center  of  mass  plus  the  rotational  motion  of  the  object  about  its 
center  of  mass.  The  slant-range,  r(t),  is  of  paramount  importance  in  the  translational 
motion  description,  whereas  angular  displacement,  0(f)  ,  is  vitally  important  in  the 
rotational  motion  description. 
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Figure  1.  Inverse  synthetic  aperture  radar  geometry. 


3.0  FORMATION  OF  THE  UNPROCESSED  ISAR 
SPATIAL  FREQUENCY  SPECTRUM 


Presented  in  this  section  is  the  underlying  theory  describing  the  formation  of  the 
unprocessed  ISAR  spatial  frequency  spectrum.  This  spectrum  results  from  the  detection 
of  the  radar  return  from  the  moving  aircraft  illuminated  by  a  linear  stepped-frequency 
modulated  radar  signal. 

The  moving  aircraft  is  actively  illuminated  by  the  ground-based  radar  transmitter. 
The  transmitted  radar  signal  is  a  stepped-frequency  modulated  waveform  as  depicted  in 
figure  2.  The  quantities  /  and  t  represent  temporal  frequency  and  time,  respectively. 
This  waveform  consists  of  N  bursts;  each  burst  contains  M  pulses.  Each  transmitted 
pulse  may  be  thought  of  as  a  quasi-monochromatic  spherical  wave  that  briefly 
illuminates  the  entire  aircraft.  This  gives  rise  to  a  coherent  reflected  (scattered)  return. 
The  return  pulse  is  quadrature-detected  at  the  radar  receiver  and  heterodyned  to  yield  a 
single  complex  number.  The  MN  complex  numbers  obtained  from  the  entire 
transmitted  radar  signal  may  be  organized  into  a  two-dimensional  (M  by  N  matrix) 
array  of  numbers.  This  array  represents  the  unprocessed  ISAR  spatial  frequency 
spectrum  of  the  aircraft.  To  understand  how  this  spatial  frequency  spectrum  is  generated 
from  a  quantitative  point  of  view,  we  will  next  show  analytically  how  a  single  complex 
spatial  frequency  spectral  value  arises  from  a  single  radar  pulse. 


3.1  SINGLE-PULSE  ILLUMINATION 


A  single  transmitted  radar  pulse  may  be  represented  by  a  quasi-monochromatic 
spherically  expanding  wave.  The  analytic  wave  function  representation  of  this  pulse  is 
given  by 


VT(r’>  0  =  A 


exp  [- ilnft ]  exp  [  i2n  {f/c)  r'] 


(1) 


where  A  is  the  amplitude,  /  the  temporal  frequency,  t  the  time,  f  the  distance  from 
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the  center  of  the  expanding  spherical  wave  (see  figure  1)  to  a  point  in  space,  and  c  the 
speed  of  light  in  vacuum.  The  aircraft  is  made  up  of  a  distribution  of  point  scatterers. 
Fgure  1  depicts  the  position  of  one  such  point  scatterer  located  at  the  coordinates  (x  ,  y). 
The  infinitesimal  reflected  return  from  this  point  scatterer,  as  seen  from  the  radar,  is 
given  by 


d\fR  (r\  t)  =  A 


exp  [- ilnft ]  exp  [j4tc  ( f/c )  r'] 

(n2 


p  (x',y')  dx'dy’. 


(2) 


The  function  p(x',y')  represents  the  reflectivity  density  of  the  aircraft  and  dx'dy'  is  an 
infinitesimal  differential  area  element  centered  about  the  point  ix' ,  y').  With  the  use  of  the 
principle  of  superposition,  the  radar  return  for  a  single  pulse  due  to  the  entire  aircraft  is 
obtained  by  integrating  equation  2.  We  obtain 


VR(0  =  exp  [-ilnft] 


oo  oo 

fj* 

-00—00 


(x\y’)  exp[i4n(f/c)r') 


dx'dy', 


(3) 


where  by  definition 

h  (x\  y')  s - -p(x',y').  (4) 

( r ') 

It  is  worth  mentioning  at  this  point  that  the  formation  of  the  final  focused  ISAR  image  is 
equivalent  to  estimating  the  reflectivity  density  function  p  (x',y')  or  equivalently 
h(x',y')  .  The  radar  pulse  return  is  detected  and  heterodyned  at  the  radar  receiver, 
yielding  a  single  complex  spatial  frequency  spectral  value 


OO  oo 

u  =  |  jh{x',y') 

-00-00 


exp  [  i4n  (//c)r'J 


dx'dy'. 


(5) 


For  this  analysis,  we  will  assume  that  the  aircraft  is  located  in  the  far  field  of  the  radar 
and  that  the  dimensions  of  the  aircraft  are  small  in  comparison  to  the  distance  r(t).  Using 
the  far-field  approximation  (Goodman,  1968),  we  may  approximate  the  quantity  r 
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appearing  in  the  exponent  of  equation  5  by  the  expression 


r' ~r  +  x'sinQ  +  y'cosQ.  (6) 

The  quantity  r’  appearing  in  the  denominator  of  equation  4  may  be  approximated  by 
r  when  using  the  far-field  approximation.  Equation  5  can  be  rewritten  as 


U  =  exp[i4n(f/c)r ]  H, 


where  by  definition 


oo  oo 


—00—00 


U',y') 


exp[i2K(px'  +  qy')} 


dx'dy' 


and 


(7) 


(8) 


p  =  2(f/c)sinQ  and  q  =  2  (f/c)  cosQ.  (9) 

It  is  important  to  note  that  the  quantity  U  appearing  in  equation  7  represents  a  single 
spatial  frequency  spectral  value  only.  In  this  analysis,  the  aircraft  may  be  regarded  as 
being  stationary  during  the  time  interval  in  which  the  aircraft  is  illuminated  by  the 
single  pulse.  Hence,  the  quantities  r(t)  and  0(0  are  effectively  constant  in  value 
over  this  time  interval. 

3.2  STEPPED-FREQUENCY  RADAR  ILLUMINATION 

The  preceding  results  will  now  be  used  to  derive  a  mathematical  expression  for  the 
two-dimensional  spatial  frequency  spectrum  generated  by  a  stepped-frequency  radar 
signal.  Keep  in  mind  that  between  radar  pulses  the  aircraft  has  moved  to  a  new  position, 
hence  r(t)  and  0(0  take  on  new  values  for  each  pulse.  During  the  time  interval 
required  to  generate  one  complete  frame  (M  by  N  array)  of  spatial  frequency  spectral 
data,  the  slant-range,  r(t),  and  the  angular  displacement,  0(0/  may  be  approximated 
by  the  first  few  terms  of  their  MacLaurin  series  expansions.  In  particular. 


r(0  =ro  +  zM+  (l/2)flof2+  (1/6) 

(10a) 

0(f)  =0  +0)  t  +  (1/2) a  f2. 

0  0  0 

(10b) 
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The  time  t  =  0  represents  the  beginning  of  the  frame.  The  translational  kinematic 
quantities  rQ,  v  ,  aQ,  and  jQ  represent  the  initial  values  of  the  slant-range,  velocity, 
acceleration,  and  jerk  respectively,  of  the  center  of  mass  of  the  aircraft  at  time  t  =  0. 
Likewise,  the  rotational  kinematic  quantities  0^,  co^,  and  represent  the  initial  values 
of  the  angular  displacement,  velocity,  and  acceleration,  respectively,  of  the  aircraft  about 
its  center  of  mass  at  t  =  0.  Let  the  quantity  t  represent  the  time  the  aircraft  is 
illuminated  by  the  mth  pulse  of  the  nth  burst  [henceforth  denoted  (m,n)th  pulse], 
where 


tmn  =  [  (m-1)  +  (n-l)M]  At. 


(ID 


The  indices  m  and  n  take  on  the  integer  values  m  =  2,  2, 3,...,  M  and  n  =  l,2,3,..., 
N.  The  constant  At  represents  the  time  interval  between  adjacent  pulses.  Substitution 
of  equation  11  into  equations  10a  and  10b  gives  the  following  expressions  for  the 
translational  and  rotational  displacements  of  the  aircraft  at  the  time  it  is  illuminated  by 
the  (m,n)th  pulse 


1  t 

r  ~  r  +v  t  +  -  a  r 
mn  o  o  mn  2  o  mn 


+  6  icfmn 


(12a) 


0  *  0  +oj  t  +^af2. 

mn  o  o  mn  2  o  mn 


(12b) 


With  the  use  of  equation  7,  the  (m,n)th  pulse  gives  rise  to  a  single  complex  spatial 
frequency  spectral  value  given  by 


U(m,n)  =  exp[i4n(fm/c)rmn)  H(m,n), 
where  with  the  use  of  equations  8  and  9 


(13) 


H  (m,  n) 


oo  oo 

{*(*•,  y> 

—00—00 


exp[i2n(pmnx'  +  qmny')]  dx'dy' 


(14) 


8 


and 


p  =2(f  /c)sinQ  and 

rmn  3m  mn 


a  =2  (/  /c)  cos0 
nmn  !m  '  mn 


(15) 


The  frequency  of  the  (m,n)th  pulse  is  given  by 


/m  =  /0+(m-DA/. 


(16) 


The  quantities  f  and  A /  represent  the  base  and  step  frequencies  of  the  linear 
stepped-frequency  transmitting  radar.  Hereafter,  H(m,n )  will  be  referred  to  as  the  ideal 
translational  motion-compensated  ISAR  spatial  frequency  spectrum  and  U(m,n)  will 
be  referred  to  as  the  unprocessed  ISAR  spatial  frequency  spectrum. 

4.0  TRANSLATIONAL  MOTION  COMPENSATION 

(TMC)  ALGORITHM 


Presented  in  this  section  is  a  description  of  the  translational  motion  compensation 
(TMC)  algorithm.  The  objective  of  the  TMC  algorithm  is  to  form  an  estimate,  H(m,n) , 
of  the  ideal  translational  motion-compensated  ISAR  spatial  frequency  spectrum,  H(m,  n), 
from  the  unprocessed  ISAR  spatial  frequency  spectrum,  U  (m,  n) .  To  do  so  requires 
estimating  the  translational  kinematic  quantities  vQ,  aQ,  and  j  . 

Equation  13  represents  the  desired  mathematical  expression  for  the  (m,n)th  element 
of  the  unprocessed  ISAR  spatial  frequency  spectrum.  This  observable  spectrum,  in  any 
real  application,  serves  as  the  input  data  to  the  ISAR  processing  system.  Substitution  of 
equation  12a  into  equation  13  gives 


U  (m,  n)  =  exp 


1  7 

i'4k  (/  /c)  (r  +v  t  +  -a  r 
x,m  o  o  mn  2  o  mn 


1  .  3 

6  Vmh 


H  (m,  n) .  (17) 


At  this  point,  lets  indicate  which  quantities  in  equation  17  are  known  and  which 
quantities  are  to  be  determined.  The  spectral  values  U(m,n )  are  the  known  observable 
inputs  to  the  ISAR  processor.  The  stepped-frequency  radar  characteristics  are  known, 
hence  /  and  ^mn  are  also  known.  The  uncooperative  aircraft's  translational 
kinematic  quantities  r  ,  vq,  aQ,  and  j  are  unknown.  If  these  quantities  were  known. 
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then  it  would  be  a  straightforward  procedure  to  solve  equation  17  for  the  ideal  spectral 
values  H(m,n).  The  objective  of  the  TMC  algorithm  is  to  estimate  the  kinematic 
quantities  vq,  aQ,  and  j  ,  which  in  turn  allows  us  to  form  estimates  of  the  spatial 
frequency  spectral  values,  denoted  by  H(m,  n) .  The  quantities  H(m,n )  would  then 
serve  as  the  inputs  to  the  RMC  algorithm.  Again,  the  RMC  algorithm  is  used  to  further 
improve  the  image  sharpness,  yielding  a  highly  focused  final  image  i  ( x ,  y)  .  These 
ideas  are  depicted  in  the  flow  diagram  shown  in  figure  3. 

As  previously  indicated,  the  principal  objective  of  the  TMC  algorithm  is  to  estimate 
the  three  translational  kinematic  quantities  vq,  aQ,  and  jQ .  To  estimate  vq,  aQ,  and  / 
we  first  multiply  the  unprocessed  spectral  values  U(m,n)  by  a  spatial  frequency  filter 
function  S(m,  n;v,a,]) .  The  quantities  v,  a,  and  j  appearing  in  the  filter  function  are 
initial  guesses  of  the  translational  velocity,  acceleration,  and  jerk,  respectively.  This  gives 


H(m,n;v,a ,])  =  S  (m,  rr,v,  a,])  U  (m,  n) , 
where  the  spatial  frequency  filter  function  is  mathematically  given  by 


(18) 


S  (m,  n;v,a,j )  =  exp 


-i4n(f  /c)  (dt  +^af2  +i?f 3  ) 
Jm  '  K  mn  2  mn  (,J  mnJ 


(19) 


The  next  step,  currently  used  in  the  TMC  algorithm,  is  to  compute  the  two-dimensional 
discrete  Fourier  transform  (2D-DFT),  (Gonzalez  and  Wintz,  1987),  of  H(m,n ;v,  a,  j)  .  In 
actual  ISAR  computer  implementations,  the  two-dimensional  fast  Fourier  transform 
(2D-FFT)  algorithm  is  employed.  This  gives 


M-lN-l 

ft(x,y;v,a,~j)  =  (m,  n;v,  a,  j)  exp  [i2n  {x  (m/M)  +y  (n/N)  }  ] .  (20) 

m  =  On  =  0 

The  function  h(x,y;v,a ,  ])  represents  an  initial  estimate  of  the  translational  motion- 
compensated  ISAR  image.  The  quality  of  this  image  may  be  determined  using  entropy 
as  a  measure  of  sharpness.  The  entropy  of  the  image  h  (x,  y;v,  a,j )  is  computed  using 
the  equation 


M-lN-l 

E  ($,&,])  =  yj h(x,y;v,  a,])\  log\h(x,y;v,a,])\.  (21) 

x  =  0  y  =  0 
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Figure  3.  ISAR  processing  algorithms. 
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With  the  help  of  a  search  algorithm,  this  process  is  iteratively  repeated  a  number  of 
times  using  different  values  of  i),  a,  and  j  until  a  value  of  velocity,  v ^  ,  acceleration, 
ag  ,  and  jerk,  ,  are  obtained  that  give  a  minimum  value  for  the  entropy. 

The  values  v  ,  a£,  and  ;  are  then  used  as  the  best  estimates  of  the  target's  actual 
translational  velocity,  acceleration,  and  jerk.  Once  these  best  estimates  have  been 
determined,  then  the  best  intermediate  translational  motion-compensated  ISAR  image 
(in  terms  of  sharpness)  is  the  one  given  by  h  (x,  y) ,  where  by  definition 

h(x,y)=hx,ywe,ae,je).  (22) 

Also, 

H(m,n)  =H(m,rr,ve,ae,je) .  (23) 

Note  that  the  phase  term  in  equation  17  containing  the  quantity  r  does  not  affect  the 
overall  sharpness  of  the  image  h  (x,  y)  .  It  does,  however,  affect  the  overall  positioning 
(or  centering)  of  the  image  (in  terms  of  translational  shift).  Hence,  if  positioning  is  of 
importance,  then  an  estimated  value  of  r g  can  be  easily  chosen  to  accomplish  a 
repositioning. 

An  alternative  technique  for  estimating  the  velocity,  acceleration,  and  jerk  is  to  use  a 
measure  where  the  computations  can  be  performed  strictly  in  the  spatial  frequency 
domain.  This  would  have  the  advantage  of  eliminating  the  numerous  2D-FFTs,  which 
must  presently  be  performed  using  the  entropy  measure,  in  estimating  the  translational 
kinematic  quantities.  The  authors  have  been  actively  investigating  the  use  of  a  number 
of  different  measures.  Based  on  preliminary  simulations,  the  burst  derivative  measure 
shows  a  great  deal  of  promise  as  an  alternative  means  for  obtaining  these  translational 
kinematic  estimates.  The  burst  derivative  measure  is  determined  by  the  equation 


B(v,  a,j)  = 


M-lN-l 

SI 

m  =  On  =  0 


~  H  (m,  n;v,  a,  ]) 


(24) 


As  in  the  case  of  entropy,  the  burst  derivative  measure  has  been  employed  successfully 
with  an  iterative  search  algorithm  to  find  estimates  ofug,  and  /  .  Once  these 
estimates  are  obtained,  then  only  one  2D-FFT  is  required  to  form  the  translational 
motion-compensated  ISAR  image  h(x,y)  from  its  estimated  translational  motion- 
compensated  spectrum  H(m,n )  . 
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5.0  TMC  ALGORITHM  SIMULATION  RESULTS 


Preliminary  comparisons  between  the  burst  derivative  measure  and  the  entropy 
measure  in  estimating  the  translational  kinematic  quantities  of  an  airborne  target  have 
been  made.  These  comparisons  have  been  done  with  simulated  as  well  as  actual  ISAR 
data  using  both  a  Sun  4  SPARCstation  2  and  a  parallel  processing  transputer  array. 
Results  typifying  the  comparison  for  simulated  data  on  the  Sun  4  SPARCstation  2  are 
presented  in  this  section. 

Let's  begin  by  referring  to  figure  4,  which  typifies  only  one  of  many  simulations 
performed.  Figure  4a  is  a  gray-scale  plot  of  the  simulated  two-dimensional  reflectivity 
density  function  p  (x',y')  of  a  fictitious  aircraft.  The  dimensions  of  this  plot 
correspond  to  76  by  76  meters.  Each  of  the  point  scatterers  depicted  in  this  figure  have 
equal  reflectivities.  For  this  simulation,  the  translational  kinematic  quantities  (ground 
truth)  describing  the  motion  of  the  aircraft  at  the  beginning  of  the  frame  were  given  by 

rQ  -  4400  m  vq  =  -4.02  m/s  a q  =  +14.20  m/s2  j  =  +0.05  m/s3, 

and  the  rotational  kinematic  quantities  (ground  truth)  were  given  by 

0^  =  90.91  deg  a>o  =  -  3.25  deg/s  a0  =  ~  0-01  deg/s2. 

The  pertinent  radar  parameters  used  for  this  simulation  were  given  by 

M  =  64  N  =  64  At  =  1.59  ms  /  =  0.372  GHz  Af  =  2.130  MHz. 

Figure  4b  is  the  gray-scale  plot  of  the  unprocessed  ISAR  spatial  frequency  spectrum 
(modulus  only  plotted).  Figure  4c  is  the  gray-scale  plot  of  the  unprocessed  ISAR  image 
one  would  obtain  by  taking  the  2D-FFT  of  the  unprocessed  ISAR  spatial  frequency 
spectrum.  If  the  unprocessed  ISAR  spatial  frequency  spectrum  is  multiplied  by  the 
spatial  frequency  filter  given  by  equation  19  for  v  =  vQ,d  =  0,;  =  0  ,  then  the  ISAR 
image  that  one  would  obtain,  from  the  2D-FFT  operation  on  the  resulting  filtered 
spectrum  is  shown  in  figure  4d.  Likewise,  if  the  unprocessed  ISAR  spectrum  is 
multiplied  by  the  spatial  frequency  filter  for  f>  =  v  ,  d  =  a  ,j  =  0,  then  the  ISAR 
image  one  would  obtain  from  the  2D-FFT  operation  on  the  resulting  filtered  spectrum  is 
shown  in  figure  4e.  Notice  in  figure  4e  that  the  aircraft  now  is  readily  discernible.  Finally, 
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Polar  Reformatting 


(g)  224 

Figure  4.  Motion  compensation  of  simulated  data  using  ground  truth. 
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if  the  unprocessed  ISAR  spectrum  is  multiplied  by  the  spatial  frequency  filter  for 
t  =  vq,  &  =  aQ,j  =  /  ,  then  the  ideal  translational  motion-compensated  ISAR  image, 
h(x,y),  one  would  obtain  from  the  2D-FFT  operation  on  the  resulting  filtered  spectrum  is 
shown  in  figure  4f .  Figure  4f  represents  about  as  well  as  one  could  hope  for  in  processing 
the  unprocessed  ISAR  spectrum  using  a  spatial  frequency  filter  containing  the  actual 
ground-truth  translational  kinematics  of  the  moving  aircraft. 

For  comparative  purposes,  if  the  translational  motion-compensated  image  obtained 
in  figure  4f  is  polar-reformatted  using  the  RMC  algorithm  employing  the  ground  truth 
rotational  kinematic  quantities  listed  on  the  previous  page,  then  the  final  ideal  rotational 
motion-compensated  (focused)  ISAR  image,  i(x,  y),  obtained  is  shown  in  figure  4g.  The 
results  shown  in  figure  4  will  serve  as  a  basis  of  comparison  for  similar  results  obtained 
using  the  entropy  measure  and  the  burst  derivative  measure  in  both  the  TMC  and  RMC 
algorithms  for  estimating  the  kinematic  quantities.  Listed  below  each  subplot  in  figure  4 
are  the  corresponding  entropy  values  for  that  subplot.  Note,  the  smaller  the  value  of  the 
entropy,  the  sharper  the  image. 

Let's  now  consider  the  use  of  the  entropy  measure  in  the  TMC  algorithm  for 
estimating  the  translational  kinematic  quantities  of  the  aircraft.  Figure  5  shows  velocity, 
acceleration,  and  jerk  slices  through  the  entropy  surface  described  by  the  entropy 
function  E(vfl,))  defined  in  equation  21.  In  particular,  figure  5a  represents  E  (v,  0, 0) 
plotted  against  v,  figure  5b  represents  E  (v£,  a,  0)  plotted  against  a,  and  figure  5c 
represents  E(t7g,  ag,/)  plotted  against/.  The  quantities  v£,  a£,  and  j  represent  the 
estimates  of  the  velocity,  acceleration,  and  jerk  using  the  entropy  measure.  These 
estimates  correspond  to  the  minima  of  figures  5a,  5b,  and  5c,  respectively.  The  triple  x 
symbol  in  these  subplots  represents  the  position  of  the  ground  truth  values.  Those 
values  of  velocity,  acceleration,  and  jerk  that  gave  rise  to  the  minimum  value  of  entropy 
were  v  =  -4.1  m/s  ,  a  =  14.4  m/s 2 ,  and  j  =  0.00  m/s3. 

Now  consider  the  use  of  the  burst  derivative  measure  in  the  TMC  algorithm  for 
estimating  the  translational  kinematic  quantities  of  the  aircraft.  Figure  5  also  shows 
velocity,  acceleration,  and  jerk  slices  through  the  burst  derivative  surface  described  by 
the  burst  derivative  function  B(vA,j)  defined  in  equation  24.  Figure  5d  represents 
6(17,0,0)  plotted  against  v,  figure  5e  represents  B(ve,a,  0)  plotted  against  a,  and 
figure  5f  represents  6  (i;^,  /)  plotted  against  j.  The  quantities  vg,  a£,  and  /  in  this 

case  represent  the  estimates  of  the  velocity,  acceleration,  and  jerk  using  the  burst 
derivative  measure.  These  estimates  correspond  to  the  minima  of  figures  5d,  5p,  and  5f, 
respectively.  Again,  the  triple  x  symbol  represents  the  position  of  the  ground  truth 
values.  Those  values  of  velocity,  acceleration,  and  jerk  that  gave  rise  to  the  minimum 
value  of  burst  derivative  were  vg  =  -3.8  m/s  ,  a e  =  14.2  m/s2  ,and  /  =  0.01  m/s3. 
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Figure  5.  Translational  velocity,  acceleration,  and  jerk  surface  slices. 


Figure  6  shows  the  resulting  translational  motion-compensated  ISAR  images 
obtained  using  both  the  entropy  and  burst  derivative  measure  kinematic  estimates  in  the 
spatial  frequency  filter  function  described  by  equation  19.  Figure  6a  shows  once  again 
the  aircraft  reflectivity  density  function  used  in  this  simulation.  Figure  6b  corresponds  to 
figure  4f  and  again  represents  the  ideal  translational  motion-compensated  ISAR  image 
obtained  using  the  ground  truth  values  in  the  spatial  frequency  filter.  Figures  4c  and  4d 
represent  the  translational  motion-compensated  ISAR  images  obtained  using  the 
entropy  and  burst  derivative  measure  kinematic  estimates,  respectively,  in  the  spatial 
frequency  filter.  Listed  below  each  subplot  in  figure  6  are  the  corresponding  entropy 
values  for  that  subplot.  We  see  that  the  ISAR  images  obtained  using  the  entropy  measure 
and  the  burst  derivative  measure  are  comparable  in  quality  and  sharpness  to  the  one 
obtained  using  the  ground  truth  information. 

In  an  actual  ISAR  imaging  application,  many  contiguous  (perhaps  overlapping) 
frames  of  data  of  a  moving  aircraft  must  be  processed.  From  practical  experience  we 
have  found  that  the  entropy  measure  must  be  employed  in  the  TMC  algorithm  to 
initially  process  the  first  frame  of  data.  This  is  due  to  the  fact  that  the  entropy  surface  has 
a  unique  global  minimum.  The  burst  derivative  surface,  on  the  other  hand,  is  periodic 
with  velocity  and,  therefore,  does  not  have  a  unique  global  minimum.  However,  once 
the  translational  kinematics  quantities  for  the  first  frame  have  been  estimated  using  the 
entropy  measure,  then  the  burst  derivative  measure  can  be  used  subsequently  in  the 
TMC  algorithm  to  track  these  kinematic  quantities  and  update  them  as  we  process 
additional  frames  of  data  in  the  sequence. 

The  TMC  algorithm  has  been  programmed  on  a  parallel  transputer  array  processor 
as  well.  One  very  important  result  of  that  work  is  that  for  64-  by  64-pixel  array  sizes  we 
have  found  that  the  burst  derivative  measure  yields  translational  motion-compensated 
images  of  comparable  image  quality  to  those  obtained  using  the  entropy  measure,  but 
does  so  in  approximately  one-sixth  the  time. 
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Figure  6.  Comparison  of  translational  motion-compensated  ISAR  images. 


6.0  ROTATIONAL  MOTION  COMPENSATION 
(RMC)  ALGORITHM 


Presented  in  this  section  is  a  description  of  the  rotational  motion  compensation 
(RMC)  algorithm.  The  objective  of  the  RMC  algorithm  is  to  compute  an 
estimate,  i  (x,  y) ,  of  the  ideal  final  focused  ISAR  image,  i(x,y),  by  using  the  output  from 
the  TMC  algorithm,  namely  the  estimated  translational  motion-compensated  ISAR 
spatial  frequency  spectrum  H(m,n)  .  To  do  so  requires  estimating  the  rotational 
kinematic  quantities  co^and  .  Please  review  figure  3  once  again. 

The  RMC  algorithm  employs  the  method  of  polar  reformatting  to  achieve  this 
objective.  In  addition,  the  RMC  algorithm  in  its  present  form  relies  heavily  on  the  use  of 
numerous  2D-FFTs  needed  in  the  entropy  measure  calculation  for  determining  the  image 
sharpness.  Preliminary  results  of  this  investigation  indicate  that  the  burst  derivative 
measure  employed  in  the  TMC  algorithm  may  be  used  in  lieu  of  the  entropy  measure  in 
the  RMC  algorithm  as  well.  This  would  have  the  advantage  of  eliminating  the  numerous 
2D-FFTs  presently  required  in  the  entropy  measure  calculations. 

To  understand  the  essential  workings  of  the  RMC  algorithm,  we  must  first  briefly 
discuss  the  essential  features  of  polar  reformatting.  The  concept  of  polar  reformatting  is 
described  in  detail  by  both  Ausherman  et  al.  (1984)  and  Wehner  (1987).  Let's  initially 
describe  the  situation  for  the  case  when  the  ideal  translational  motion-compensated 
ISAR  spatial  frequency  spectrum,  H(m,n)  ,  is  known. 

Recall  that  the  2D-FFT  of  H(m,n)  led  to  the  ideal  translational  motion-compensated 
ISAR  image,  h(x,y),  shown  in  figure  4f.  By  polar  reformatting  the  spectral  values  H(m,n), 
the  ideal  final  focused  ISAR  image,  i(x,y),  shown  in  figure  4g  was  obtained.  Let's  now 
describe  the  necessary  steps  leading  to  the  polar-reformatted  image  shown  in  figure  4g. 

For  the  two-dimensional  scenario  being  discussed,  the  spectral  values  H  (m,n)  are 
physically  located  on  a  two-dimensional  spatial  frequency  plane.  The  Cartesian 
coordinates  describing  the  position  of  each  value  of  H(m,n )  are  specified  by  the 
ordered  pair  of  numbers  (pmn,qmn)  as  determined  by  equation  15.  These  spectral 
value  coordinates  in  turn  depend  on  both  fm  and  0  as  indicated  in  equations  16  and 
12b,  respectively.  The  angles  0  depend  on  the  ground  truth  rotational  kinematic 
quantities  0  >  G>  >  and  .  Using  the  ground  truth  rotational  kinematic  values  and 
radar  parameters  listed  on  page  13  leads  to  the  two-dimensional  spatial  frequency 
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domain  plots  shown  in  figure  7.  Figure  7a  depicts  a  Cartesian  set  of  axes  labelled  p  and  q. 
Each  value  of  (p  ,  9mn)  obtained  from  equation  15  is  plotted  as  a  point  (black  dot)  in 
this  plane.  The  wedge-shaped  grid  pattern  resulting  in  figure  7a  represents  the  region 
where  the  64-  by-64  array  of  ideal  spectral  values  are  physical  located  for  the  simulation 
considered  in  this  technical  report.  A  magnified  view  of  this  wedge-shaped  grid  pattern 
is  shown  in  figure  7b,  where  the  individual  black  dots  comprising  the  array  are  now 
easily  distinguishable.  In  using  the  2D-FFT  to  construct  the  translational  motion- 
compensated  image  in  figure  4f  from  H(m,n),  it  was  tacitly  assumed  that  the  spectral 
values  H(m,n)  were  located  on  a  rectangular  grid  pattern,  and  not  on  a  wedge-shaped 
grid  pattern  as  seen  in  figure  7. 

Polar  reformatting  is  a  technique  used  to  construct  a  new  array  of  spectral  values 
positioned  on  a  rectangular  grid  pattern  from  an  array  of  spectral  values  lying  on  a 
wedge-shaped  grid  pattern.  This  new  array  of  spectral  values  when  Fourier  transformed 
will,  in  general,  yield  a  sharper  (more  highly  focused)  ISAR  image  in  comparison  to  the 
one  obtained  when  polar  reformatting  is  not  used.  The  polar  reformatting  technique 
used  for  the  simulations  performed  in  the  Processing  Research  and  Development  Branch 
(Code  761)  of  the  NCCOSC  RDT&E  Division  involves  a  five-step  approach: 

1.  Oversample  the  64-by-64  array  of  spectral  values  by  a  factor  of  four.  This  yields  a 
new  256-by-256  array  of  spectral  values. 

2.  Rotate  the  newly  formed  wedge-shaped  pattern  of  256  by  256  points  obtained  in 
step  1  until  the  wedge  bisector  lies  along  the  positive  p-axis;  see  figure  8a. 

3.  Overlay  a  64-by-64  rectangular  grid  pattern  of  points  of  maximum  extent,  but 
wholly  contained  within  the  newly  formed  wedge-shaped  grid  pattern  obtained  in 
step  2;  see  figure  8b. 

4.  At  each  point  on  the  rectangular  grid  pattern  generated  in  step  3,  find  the  three 
nearest  neighboring  points  on  the  wedge-shaped  grid  pattern. 

5.  At  each  point  in  the  rectangular  array,  determine  a  new  spectral  value  by 
interpolating  the  spectral  values  associated  with  the  three  nearest  neighbor  points 
found  in  step  4. 

Applying  this  five-step  process  gives  rise  to  a  64-by-64  rectangular  array  of  spectral 
values  whose  discrete  Fourier  transform  yields  the  desired  focused  ISAR  image,  i(x,  y), 
shown  in  figure  4g. 
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Figure  7.  Position  of  the  ideal  spatial  frequency  spectral  values. 
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The  ideal  translational  motion-compensated  spectral  values,  H{m,n)  ,  are,  of 
course,  unknown.  Hence,  the  estimated  spectral  values,  obtained  from  the 

TMC  algorithm  must  be  used  instead.  In  addition,  the  ground  truth  rotational  kinematic 
quantities  of  the  moving  aircraft  are  also  unknown  and  must  be  estimated.  The 
following  procedure  summarizes  the  key  steps  of  the  RMC  algorithm. 

1.  Guess  initial  values,  (0  and  a  ,  of  the  ground  truth  angular  velocity,  to  ,  and  the 
angular  acceleration,  a^,  respectively.  Let  the  estimate,  0 ,  of  the  ground  truth 
rotation  angle,  0^,  be  zero. 

2.  Use  the  values  from  step  1  with  the  known  values  of  t  to  generate  the  quantities 

®mn  using  equation  12b. 

3.  Use  the  values  from  step  2  and  the  quantities,  /  ,  from  equation  16  to  generate  the 

ordered  pairs  using  equation  15. 

4.  Using  the  order  pairs  from  step  3,  form  a  wedge-shaped  grid  pattern  analogous  to 
the  one  formed  in  figure  7  using  actual  ground  truth  data. 

5.  Apply  the  five-step  reformatting  process  described  on  the  preceding  page  to  the 
wedge-shaped  grid  pattern  obtained  in  step  4  above.  This  gives  an  estimated 
rotational  motion-compensated  ISAR  spatial  frequency  spectrum,  denoted  by 

2(m,  n;cb,  a)  . 

6.  Perform  a  2D-FFT  on  the  result  obtained  in  step  5.  This  gives  an  estimated 
rotational  motion-compensated  ISAR  image,  denoted  by  i(x,  y;co,  a)  . 

7.  Calculate  the  entropy  of  the  image  obtained  from  step  6  using  the  following 
equation: 

M-lN-l 

£(<o,a)  =  y  yji(s,y ;<&,&)[  /og|i  (x,  y;<6,  a)|.  (25) 

*  =  0  y  =  0 

At  this  point,  new  estimates  of  the  rotational  kinematic  quantities  are  chosen  in  step  1, 
and  the  above  procedure  is  repeated.  This  process  is  repeated  over  and  over  again,  using 
a  search  algorithm  similar  to  the  one  used  in  the  TMC  algorithm,  until  a  minimum  value 
for  the  entropy  is  obtained.  The  estimated  rotational  kinematic  quantities,  co^  and  a  , 
giving  the  minimum  value  of  entropy  are  used  as  the  best  estimates  of  the  ground  truth 
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angular  velocity  and  acceleration.  The  above  procedure  is  applied  one  last  time  using  the 
rotational  kinematic  estimates,  <og  and  ,  to  obtain  the  best  estimate,  i  (x,  y)  ,  of 
the  ideal  final  focused  ISAR  image  i(x,y),  where  by  definition 


i(x,y)  =t(x,y;  ©g>ap. 

Also, 

2  (m,  n )  =2  (m,  «;cog,  ap . 

By  using  the  burst  derivative  measure  (see  equation  28  below),  we  can  replace  steps  6 
and  7  in  the  RMC  algorithm  by  calculating  the  burst  derivative  measure  of  the  polar- 
reformatted  spectrum  obtained  in  step  5  and  thus  eliminate  the  multitude  of  2D-FFTs 
presently  used  when  the  entropy  measure  is  used  in  measuring  image  sharpness. 


(26) 


(27) 


B  (<&,  a) 


M-lN-l 

m  =  0n  =  0 


(28) 


As  in  the  case  of  entropy,  the  burst  derivative  measure  has  been  employed  successfully 
with  an  iterative  search  algorithm  to  find  estimates  of  co^  and  .  Once  these 
estimates  are  obtained,  then  only  one  2D-FFT  is  required  to  form  the  rotational  motion- 
compensated  ISAR  image  i  (*,  y)  from  its  estimated  rotational  motion-compensated 
spectrum  2(m,  n)  . 

7.0  RMC  ALGORITHM  SIMULATION  RESULTS 

Preliminary  comparisons  between  the  burst  derivative  measure  and  the  entropy 
measure  in  estimating  the  rotational  kinematic  quantities  of  an  airborne  target  have  been 
made.  These  comparisons  have  been  done  with  simulated  as  well  as  actual  ISAR  data 
using  a  Sun  4  SPARCstation  2.  Results  typifying  the  comparison  for  simulated  data  on 
the  Sun  4  SPARCstation  2  are  presented  in  this  section.  In  particular,  the  translational 
and  rotational  kinematic  quantities  (ground  truth)  as  well  as  radar  parameters  used  in 
the  simulation  described  on  page  13  in  section  5  apply  here  as  well.  For  convenience. 
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these  values  are  listed  below: 


r  =  4400  m  v  =  -4.02  m/s  a  =  +14.20  m/s2  i  =  +  0.05  m/s3 

O  O  O  1 0 

=  90.91  deg  (0^  =  -  3.25  deg/s  a0  =  ~  deg/s2 
M  =  64  N  =  64  At  =  1.59  ms  f  =  0.372  GHz  A f  =  2.130  MHz. 

First  consider  the  use  of  the  entropy  measure  in  the  RMC  algorithm  for  estimating 
the  rotational  kinematic  quantities  of  the  aircraft.  Figure  9  shows  angular  velocity  and 
angular  acceleration  slices  through  the  entropy  surface  described  by  the  entropy 
function  E  (co,  a)  defined  in  equation  25.  In  particular,  figure  9a  represents  E  (to,  0) 
plotted  against  to,  and  figure  9b  represents  E  (co^,  a,  0)  plotted  against  a.  The 
quantities  (o^  and  ag  represent  the  estimates  of  the  angular  velocity  and  acceleration 
using  the  entropy  measure.  These  estimates  correspond  to  the  minima  of  figures  9a  and 
9b,  respectively.  The  triple  x  symbol  in  these  subplots  represents  the  position  of  the 
ground  truth  values.  Those  values  of  angular  velocity  and  acceleration  that  gave  rise  to 
the  minimum  value  of  entropy  were  0)  =-3.20  deg/s  and  a  =  0.00  deg/s 1 . 

Now  consider  the  use  of  the  burst  derivative  measure  in  the  RMC  algorithm  for 
estimating  the  rotational  kinematic  quantities  of  the  aircraft.  Figure  9  also  shows  angular 
velocity  and  acceleration  slices  through  the  burst  derivative  surface  described  by  the 
burst  derivative  function  B  (co,  a)  defined  in  equation  28.  Figure  9c  represents 
B(o),  0)  plotted  against  co,  and  figure  9d  represents  B(o)^,  a)  plotted  against  a.  The 
quantities  co^  and  ag  in  this  case  represent  the  estimates  of  the  angular  velocity  and 
acceleration  using  the  burst  derivative  measure.  These  estimates  correspond  to  the 
minima  of  figures  9c  and  9d,  respectively.  Again,  the  triple  x  symbol  represents  the 
position  of  the  ground  truth  values.  Those  values  of  angular  velocity  and  acceleration 
that  gave  rise  to  the  minimum  value  of  burst  derivative  were  co  =  -  2.30  deg/s  and 
ag  =  -  0.01  deg/s2  . 

Figure  10  shows  the  resulting  rotational  motion-compensated  ISAR  images  obtained 
using  the  entropy  and  burst  derivative  measure  kinematic  estimates.  Figure  10a  shows 
once  again  the  aircraft  reflectivity  density  function  used  in  this  simulation.  Figure  10b 
corresponds  to  figure  4g  and  again  represents  the  ideal  rotational  motion-compensated 
ISAR  image  obtained  using  the  ground  truth  values.  Figures  10c  and  lOd  represent  the 
rotational  motion-compensated  ISAR  images  obtained  using  the  entropy  and  burst 
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ENTROPY 


BURST  DERIVATIVE 


Angular  Velocity  Slice  Angular  Velocity  Slice 


-8.0  ©  -1.0  -8.0  to  -1.0 


Angular  Acceleration  Slice  Angular  Acceleration  Slice 


-  0.80  a  +  0.80  -  0.80  a  +  0.80 


Figure  9.  Angular  velocity  and  acceleration  surface  slices. 
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Reflectivity  Function 


(a) 

Ideal  Ground  Image 


(b)  224 

Burst  Derivative  Processed  Image 


(d)  239 


Entropy  Processed  Image 


(c)  221 


Figure  10.  Comparison  of  rotational  motion-compensated  ISAR  images. 


27 


derivative  measure  kinematic  estimates,  respectively.  Listed  below  each  subplot  in  figure 
10  are  the  corresponding  entropy  values  for  that  subplot.  We  see,  once  again,  that  the 
ISAR  images  obtained  using  the  burst  derivative  measure  and  the  entropy  measure  are 
comparable  in  quality  and  sharpness  to  the  one  obtained  using  the  ground  truth 
information. 

8.0  TMC  AND  RMC  RESULTS  FOR  A  REAL  AIRCRAFT 

The  results  presented  in  sections  5.0  and  7.0,  involving  the  TMC  and  RMC 
algorithms,  respectively,  are  typical  of  simulations  performed  to  date  on  a  Sun  4 
SPARCstation  2  in  the  Processing  Research  and  Development  Branch  (Code  761)  to  test 
the  validity  of  the  burst  derivative  measure  in  its  ability  to  predict  both  the  translational 
and  rotational  kinematics  of  an  uncooperative  airborne  target.  In  this  section,  we  will 
present  results  on  the  use  of  the  TMC  and  RMC  algorithms,  using  the  burst  derivative 
measure,  for  processing  actual  ISAR  data  of  a  DC-9  commercial  airline  flying  out  of 
Lindbergh  Field  in  San  Diego,  California.  The  ISAR  data  was  furnished  by  the  Radar 
Branch  (Code  755)  of  the  NCCOSC  RDT&E  Division.  The  pertinent  radar  parameters 
used  in  collecting  the  ISAR  data  are 


M  =  64  N  =  64  At  =  0.146  ms  f  =  9.00  GHz  A f  =  4.00  MHz. 

Figure  11  summarizes  the  results  for  one  frame  of  ISAR  data  processed  on  a  Sun  4 
SPARCstation  2.  Figure  11a  represents  an  array  dots  depicting  the  outline  of  a  DC-9 
aircraft.  This  subplot  has  nothing  to  do  with  the  actual  ISAR  data  processed.  This  figure 
is  merely  included  to  remind  the  reader  of  the  general  shape  of  a  DC-9  commercial 
airliner.  Figure  lib  is  a  gray-scale  plot  of  the  actual  DC-9  unprocessed  ISAR  spatial 
frequency  spectrum  (modulus  only  plotted).  The  actual  spectral  values  supplied 
consisted  of  a  64-by-64  array  of  complex  numbers.  Figure  11c  is  the  gray-scale  plot  of  the 
unprocessed  ISAR  image  one  obtains  by  taking  the  2D-FFT  of  the  unprocessed  ISAR 
spatial  frequency  spectrum.  When  the  unprocessed  ISAR  spatial  frequency  spectrum  is 
multiplied  by  the  spatial  frequency  filter  given  by  equation  19  for  v  =  vg,a  =  0,;  =  0  , 
the  ISAR  image  that  one  obtains  by  taking  the  2D-FFT  of  the  resulting  filtered  spectrum 
is  shown  in  figure  lid.  Likewise,  when  the  unprocessed  ISAR  spectrum  is  multiplied  by 
the  spatial  frequency  filter  for  v  =  v  ,a  =  a  ,  ]  =  0  ,  the  ISAR  image  that  one  obtains 
from  taking  the  2D-FFT  of  the  resulting  filtered  spectrum  is  shown  in  figure  lie.  Next  the 
unprocessed  ISAR  spectrum  is  multiplied  by  the  spatial  frequency  filter  for 
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DC-9  Profile 


Unprocessed  Spectrum  Unprocessed  Image 


Velocity  Compensation  Acceleration  Compensation  Jerk  Compensation 


Polar  Reformatting 


Image  Clipping 


Figure  11.  Motion  compensation  of  actual  DC-9  ISARdata 


t>  =  v^ii  =  =  jg.  The  translational  motion-compensated  ISAR  image  one  obtains 

by  taking  the  2D-FFT  of  the  resulting  filtered  spectrum  is  shown  in  figure  Ilf.  Those 
values  of  estimated  translational  velocity,  acceleration,  and  jerk  that  gave  rise  to  the 
minimum  value  of  burst  derivative  measure  are  given  by 

v£  =  -21.3  m/s  =  -2.4  m/s2  jg  =  0.00  m/s3. 

Figure  11  g  shows  the  final  rotational  motion-compensated  ISAR  image  obtained 
using  the  burst  derivative  measure  kinematic  estimates.  Those  values  of  the  estimated 
angular  velocity  and  angular  acceleration  that  gave  rise  to  the  minimum  value  of  the 
burst  derivative  measure  are  given  by 

(og  =  -3.10  deg/s  ag  =  0.00  deg/s2. 

To  enhance  the  contrast  of  figure  llg,  a  simple  image  clipping  operation  was 
performed.  The  result  is  shown  in  figure  llh.  Pixel  values  in  figure  llg  range  from  0  to 
255  (8-bit  quantization).  Those  values  below  175  and  above  215  were  set  to  zero.  After 
performing  this  clipping  operation,  the  remaining  non-zero  pixel  values  were 
normalized,  i.e.,  minimum  set  to  0  and  the  maximum  to  255.  Listed  below  each  subplot 
in  figure  11  are  the  corresponding  entropy  values  for  that  subplot. 


9.0  CONCLUSIONS 

Inverse  synthetic  aperture  radar  (ISAR)  is  an  imaging  technique  that  shows  great 
promise  in  classifying  airborne  targets  in  realtime  under  all  weather  conditions.  The 
success  of  classifying  targets  using  ISAR  is  predicated  upon  forming  highly  focused 
radar  images  of  these  targets.  Over  the  past  several  years,  a  large  body  of  ISAR  data  has 
been  collected  and  considerable  effort  has  been  expended  to  develop  algorithms  to  form 
high-resolution  focused  images  from  these  data.  Two  key  algorithms  required  in 
processing  of  ISAR  data  are  the  translational  motion  compensation  (TMC)  and  rotational 
motion  compensation  (RMC)  algorithms.  Vital  to  each  of  these  algorithms  is  the  ability 
to  measure  the  sharpness  of  an  image  during  the  iterative  computational  procedure 
required  to  process  and  focus  these  images.  At  the  present  time,  the  key  figure  of  merit 
used  to  measure  image  sharpness  is  entropy.  However,  a  new  figure  of  merit,  called  the 
burst  derivative  measure,  offers  a  reduction  in  the  computational  complexity  for  the 
TMC  and  RMC  algorithms. 
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Based  on  our  preliminary  investigation  of  the  burst  derivative  measure,  we  have 
reached  the  following  conclusions: 


1.  The  burst  derivative  measure  is  comparable  to  the  entropy  measure  in 
terms  of  its  ability  to  yield  highly  focused  ISAR  images. 

ISAR  images  have  been  processed  on  both  a  Sun  4  SPARCstation  2  as  well  as 
a  parallel  processing  transputer  array.  The  quality  and  sharpness  of  images 
obtained  using  the  burst  derivative  measure  are  comparable  to  those 
obtained  using  the  entropy  measure. 


2.  The  burst  derivative  measure  offers  a  speed  advantage  over  the  entropy 
measure. 

The  entropy  measure  determines  image  sharpness  in  real  space,  whereas  the 
burst  derivative  measure  determines  image  sharpness  in  reciprocal  space. 
The  spatial  frequency  filtering  operations  performed  by  the  TMC  algorithm 
and  the  polar  reformatting  and  interpolation  computations  performed  by 
the  RMC  algorithm  are  carried  out  in  reciprocal  space.  Hence,  the  numerous 
2D-FFTs  presently  required  to  go  from  reciprocal  space  to  real  space  when 
using  the  entropy  measure  can  be  avoided. 

When  tested  on  a  parallel  proce'-.ing  transputer  array,  it  was  found  that  the 
TMC  algorithm  using  the  burst  derivative  measure  had  a  six  times  speed 
advantage  over  the  TMC  algorithm  using  the  entropy  measure  for  64-by-64 
array  sizes. 


3.  Entropy  must  be  used  in  processing  the  first  frame  of  ISAR  data. 

Entropy  is  needed  to  initially  estimate  the  kinematics  of  the  airborne  target. 
Once  these  estimates  have  been  made,  then  the  burst  derivative  measure  can 
track  the  kinematics  as  they  dynamically  change  in  time. 
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